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Abstract 

We prove that the maximum norm of the deformation tensor of velocity 
gradients controls the possible breakdown of smooth (strong) solutions for 
the 3-dimensional compressible Navier-Stokes equations, which will hap- 
pen, for example, if the initial density is compactly supported [33]. More 
precisely, if a solution of the compressible Navier-Stokes equations is ini- 
tially regular and loses its regularity at some later time, then the loss of 
regularity implies the growth without bound of the deformation tensor as 
the critical time approaches. Our result is the same as Ponce's criterion 
for 3-dimensional incompressible Euler equations ( [26]). Moreover, our 
method can be generalized to the full Compressible Navier-Stokes system 
which improve the previous results. In addition, initial vacuum states are 
allowed in our cases. 
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1 Introduction 

The time evolution of the density and the velocity of a general viscous compress- 
ible barotropic fluid occupying a domain Q G R'^ is governed by the compressible 
Navier-Stokes equations 

dtp + div(pu) = 0, 

dt{pu) + div(pM ®u)- nAu - + A)V(divM) + VP(p) = 0, 

where p, u, P denotes the density, velocity and pressure respectively. The equation 
of state is given by 

P{p) = ap^ (a > 0, 7 > 1), 

H and A are the shear viscosity and the bulk viscosity coefficients respectively. 
They satisfy the following physical restrictions: 

2 

^>0,A + -^>0. (1.2) 

3 

The equations (11. ip will be studied with initial conditions: 

(p,m)(x,0) = (po, wo)(a:), (1.3) 

and three types of boundary conditions: 

1) Cauchy problem: 

n = and ( in some weak sense) p,u vanish at infinity; (1-4) 

2) Dirichlet problem: in this case, is a bounded smooth domain in R^, and 

u = Oon on; (1.5) 

3) Navier-slip boundary condition: in this case, i7 is a bounded smooth domain 
in R^, and 

u ■ n = 0, curlti x n = on dfl (1-6) 

where n = (ni, 7^2,^3) is the unit outward normal to dQ. The first condition 
in (II. 6p is the non-penetration boundary condition, while the second one is also 
known in the form 

{T>{u) ■ n)r = —KrUr, (1.7) 

where ^^(u) is the deformation tensor: 

V{u) = ^(yu + Vu'), (1.8) 



2 



and Hr is the corresponding principal curvature of dQ. Condition (11 .Tp implies 
the tangential component of T>{u) ■ n vanishes on flat portions of the boundary 
dfl. Note that Vn can be decomposed as 



(1.9) 



where T>{u) is the deformation tensor defined by ( 11. Sp and 




(1.10) 



known as the rigid body rotation tensor. The tensors 'C'(m) and S{u) are respec- 
tively the symmetric and skew-symmetric parts of Vu. 

There are huge literatures on the large time existence and behavior of solutions 
to (II. ip . The one-dimensional problem has been studied extensively by many peo- 
ple, see P^ [ ^ [2 ^ [5n] and the references therein. The multidimensional problem 
(II. ip was investigated by Matsumura-Nishida [21], who proved global existence 
of smooth solutions for data close to a non-vacuum equilibrium, and later by 
Hoff [HHIS] for discontinuous initial data. For the existence of solutions for arbi- 
trary data, the major breakthrough is due to Lions [22|[23] (see also Feireisl [T3]). 
where he obtains global existence of weak solutions - defined as solutions with 
finite energy - when the exponent 7 is suitably large. The main restriction on 
initial data is that the initial energy is finite, so that the density is allowed to 
vanish. 

However, the regularity and uniqueness of such weak solutions remains com- 
pletely open. It should be noted that one cannot expect too much regularity of 
Lions's weak solutions in general because of the results of Xin( [33]), who showed 
that there is no global smooth solution (p, u) to Cauchy problem for (II. ip with a 
nontrivial compactly supported initial density. Xin's blowup result ( [33]) raises 
the questions of the mechanism of blowup and structure of possible singularities: 
What kinds of singularities will form in finite time? What is the main mechanism 
for possible breakdown of smooth solutions for the 3-D compressible equations? 

We begin with the local existence of strong (or classical) solutions. In the 
absence of vacuum, the local existence and uniqueness of classical solutions are 
known [25|[3T]. In the case where the initial density need not be positive and may 
vanish in an open set, the existence and uniqueness of local strong (or classical) 
solutions are proved recently in [1HSIE1I2B]- Before stating their local existence 
results, we first give the definition of strong solutions. 



Definition 1.1 (Strong solutions) {p,u) is called a strong solution to U.l\) in 
fl X (0,T), if for some Qq G (3,6], 



o<pe C([o,T], w^i'^°(i])), pt e C([o,r],L''«(n)) 
u e c([o, T],Dl,n D\n)) n l\o, t- D^^'^^n)) 

ut G L°°(0, T,; L\^l)) f] 1^(0, T; Z^^(l^)), 



(1.11) 



and {p,u) satisfies U.l]) a.e. in Q x (0,T). 
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Here and throughout this paper, we use the following notations for the standard 
homogeneous and inhomogeneous Sobolev spaces. 

fD^^^{n) = {u e LUQ) : \\V'u\\Lrin) < oo}, 
Dl{n) = {ue L%n) : \\Vu\\L2^n) < oo, and ([HD or or holds} , 

In particular, Cho etc ^ proved the following result. 
Theorem 1.1 If the initial data po and uq satisfy 

o< poe L\n)r]W^'^{n), uoe DlnD\n), (1.12) 

for some g G (3, oo) and the compatibility condition: 

— fiAuo — (A + p)'Vdwuo + VF(po) = Po^^9 for some g G L^i^), (1-13) 

then there exists a positive time Ti G (0,oo) and a unique strong solution {p,u) 
to the initial boundary value problem ( fi.ij) ( fi.3l) together with ( |j.^[ ) or or 
l[1.6\) in Q X (0,Ti]. Furthermore, the following blow-up criterion holds: if T* is 
the maximal time of existence of the strong solution {p, u) and T* < oo, then 



sup (llpll/finwi.^o + \\u\\d},) = oo, (1.14) 

with q = min(6, q). 

There are several works ( [9l[TT | [T2 | [T7 t l20] ) trying to establish blow up criterions 
for the strong (smooth) solutions to the compressible Navier-Stokes equations. In 
particular, it is proved in [llj for two dimensions, if 7/i > 9A, then 

lim ( sup IIpIU- + / (IIpIIh-1,,0 + ||Vp||i2)rft ) = cx). 



T-5>T* 



0<t<T 



where T* < oo is the maximal time of existence of a strong solution and (/o > 3 
is a constant. 

Later, we [T71[T9],[20] first establish a blowup criterion, analogous to the Beal- 
Kato-Majda criterion [T] for the ideal incompressible flows, for the strong and 
classical solutions to the isentropic compressible flows in three-dimension: 

lim / II VwllLoorft = oo, (1-15) 

under the stringent condition on viscous coefficients: 

7p>X. (1.16) 
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Recently, the ideas of [TTllinilSQ] has been generahzed in |T2] to estabhsh a 
blowup criterion similar to fll.lSp . under the same assumption (II. 161) . for the 
non-isentropic fluids, that is, 

lim sup II 6*11 LOO + [ \\V u\\ LOO dt] = oo. (1.17) 

T^T* \o<t<T Jo J 



Very recently, in the absence of vacuum, Huang-Li [18] succeeded in removing the 
crucial condition fll.l6p of |T2[[T71[19l|20] and established blowup criterions (11.151) 
and ^ 

lim / {\\e\\l^ + \\Vu\\L^)dt = oo, (1.18) 
Jo 

for isentropic and non-isentropic compressible Navier-Stokes equations, under the 
physical restrictions (11.21) . respectively. 

It should be noted that for ideal incompressible flows, Beal-Kato-Majda [1] 
established a well-known blowup criterion for the 3-Dimensional incompressible 
Euler equations that a solution remains smooth if 

\\S{u)\\Loodt (1.19) 

is bounded, where S{u) is the rigid body rotation tensor defined by (11.101) . Later, 
Ponce p6] rephrased the Beal-Kato-Majda's theorem in terms of deformation 
tensor T>{u), that is, the same results in ^Ij hold if 

\\V{u)\\Loodt (1.20) 

remains bounded. Moreover, as pointed out by Constantin [10], the solution is 
smooth if and only if 







is bounded, where ^ is the unit vector in the direction of vorticity curlw. All 
these facts in (TUTDIESI show that the solution becomes smooth either the skew- 
symmetric or symmetric part of Vm is controlled. 

Note that the results in [IB] are not so satisfactory in two-fold: one is that 
the results exclude initial vacuum states; moreover, nothing is known from (11.151) 
about the natural question: which part of Vu, the symmetric part I'('u) or the 
skew-symmetric part S{u), will become arbitrarily large as the critical time ap- 
proaches? 

The aim of this paper is to improve all the previous blowup criterion results for 
compressible Navier-Stokes equations by removing the stringent condition (11.161) . 
and allowing initial vacuum states, and furthermore, instead of (11.151) . describing 
the blowup mechanism in terms of the deformation tensor V{u). Our main result 
can be stated as follows: 
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Theorem 1.2 Let {p, u) he a strong solution of the initial boundary value problem 
( (i. jp / (j.5j) together with or M.^) or M.(^) satisfying fli.i il) . Assume that the 
initial data {po,Uo) satisfies U.l^) and U.13\) . If T* < oo is the maximal time 
of existence, then 



where T>{u) is the deformation tensor defined by U.8\) . 
A few remarks are in order: 

Remark 1.1 Theorem also holds for classical solutions to the compressible 
flows with initial vacuum, which improves the results of IT8^ to the case where the 
initial density need not be positive and may vanish in an open set. In addition, 
Theorem M.^ holds for all p and A satisfying the physical restrictions ( (i.^j) . which 
removed the condition M.lb]) which is essential in the analysis in /7^[77| , [7P[ [^ ^ . 

Remark 1.2 In 1998, Xin fS&j gave an life span estimate of classical solutions 
to the compactly supported initial density of the Cauchy problem U.l\)U.!^)ljil.4\ ) 
at least in one dimension. However, it's unclear which quantity becomes infinite 
as the critical time approaches. Theorem \l.S\ shows that instabilities can develop 
only if the size of the deformation tensor becomes arbitrarily large. 

Remark 1.3 Theorem I i . ij| gives a counter part of Ponce's result in Wd^ for the 
incompressible flows. 

Next, we indicate that the results in Theorem 11.21 can be generahzed to the 
non-isentropic fluids described by 

(dtp + div(pM) = 0, 
dt{pu) + div(pu m) — pAu — {p + A)V(divtt) + VP = 0, q 22) 



where 6 is the absolute temperature, P = Rp9{R > 0), and > 0, i? > 0, c^, > 
are physical constants. 

The local existence of strong solutions with initial vacuum is established in [5], 
where it is essentially shown that for initial data {po,uq,6q) satisfying 




(1.21) 



Cv[dt{pO) + divipOu)] - kA9 + Pdivu 



Vu + Vu^\'^ + X{divuf 



< po e W^'^{n) for some 3 < ^ < 6 

uo e H^{n)r]H^{n),9o e H\n),eo > o, 



(1.23) 



and the compatibility conditions 



pAuo + {p + A)Vdiviio - RV{poOo) = p^gi, 

P i 

kAOo + -\Vuo + VmqP + A(divMo)^ - RpoOodivuo = p§g2 



(1.24) 
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for some gi,g2 G L'^{^), furthermore, {x G Q\po{x) = 0} being an open subset 
of Q, there exist a > and a unique strong solution {p,u,9) on [0,T^] to the 
Cauchy problem of (11.221) . such that for any go ^ (3, q), 

< p G C{[0,T,],W'''''), pt G C([0,T,],L''«), 
u G C([0,T,], n D^) n L2(0, T,; D^'-?"), 

G L~(0,r,;L2) nL2(0,T,;D^), (1.25) 
6 G C([0, T,]; i/^) n l2(0, T,; D^'^'), ^ > 0, 
0^ G L°°(0,T,;L2)nL2(0,T,;i/^). 



By modifying the analysis for Theorem 11.21 and in [TS], one can obtain the 
following blowup criterion for the full compressible Navier-Stokes system (11.221) . 

Theorem 1.3 Assume that the initial data satisfy (1 1 . 23^ and (| 1 . 24} ■ Let (p, u, 6) 



he a strong solution to the Cauchy problem of {\1.22\] satisfying {\1.25\] . If T* < oo 
is the maximal time of existence, then 

hm^ r {\\e\\l^ + \\V{u)\\l^) dt = oo. (1.26) 
Jo 

As aforementioned [3133], there are no global smooth solutions for the com- 
pressible Navier-Stokes when the initial density is compactly supported. We have 
the following corollary immediately. 

Corollary 1.4 Assume that {po,uo,6o) G H^{R'^) satisfy the initial compatibility 
condition i \1.24\> s'wc/i that there exists a finite number < r < oo with supppo C 
Br- Let {p,u,9) be the corresponding classical solutions. Then there exists a time 
< oo, such that M.26\) holds. 

We now comment on the analysis of this paper. Note that in all previous 
works [T71420] (see also [12]), either the assumption (11.161) or the absence of vacuum 
played an important role in their analysis in order to obtain an improved energy 
estimate which is essential not only for bounding the norm of the convection 
term F = put+pu-Vu but also for improving the regularity of the solutions. Their 
method depends on the L°°-norm of Vm also. It is thus difficult to adapt their 
analysis here. To proceed, some new ideas are needed. The key step in proving 
Theorem 11.21 is to derive the L^-estimate on gradients of both the density p and 
the velocity u. Observe that there are two main difficulties here: one is due to the 
possible vacuum states, the other is the strong nonlinearities of convection terms. 
In order to overcome these difficulties, we will use the simple observation that 
the momentum equations (ll.ip g become "more" diffusive near vacuum if divided 
on both sides by p as long as p remains bounded above which is guaranteed 
by the boundedness of the temporal integral of the super-norm in space of the 
deformation tensor. Thus a new energy estimate by using the effective stress 
tensor will lead to a prior estimates on the L^-norms of gradients of both the 
density and the velocity. A simple combination of the above facts with the ideas 
used in [18] then yields the blowup criterion for the full compressible Navier- 
Stokes system (ll.22p . The details of the proof of Theorems 11.21 and 11.31 are given 
in Sections 2 and 3 respectively. 
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2 Proof of Theorem 11.2 



Let (p, m) be a strong solution to the problem fll.ip - fll.2l) as described in The- 
orem [L2l First, the standard energy estimate yields 



sup {Wp'^Mml- + \\p\\h) + r W'^uWl.dt <C, 0<T<T*. (2.1) 

0<t<T Jo 



_ i\\p'^Mt)fr,+\\p\\L)+ r 

0<t<T 

To prove the theorem, we assume otherwise that 



lim [ \\V{u)\\Lo^^n)dt < C < oo. (2.2) 

T^T" Jq 

Then (12. 2p . together with (ll.ll) i . immediately yields the following L°° bound of 
the density p. Indeed, on has 

Lemma 2.1 Assume that 



[ \\divu\\L^dt <C, 0<T <T*. 
Jo 



Then 

sup IIpIIloo <C, 0<T <T*. (2.3) 

0<t<T 

Proof. It follows from (ll.ip i that for Vp > 7, 

dtipP) + divipPu) + {p- l)/divw = 0. (2.4) 
Integrating (12. 4p over Q leads to 

dt / //dx < (p - l)||divM||Loc(j^) / pPdx, 
Jn Jn 



that is. 



p — 1 

dt\\p\\LP < -^\\divu\\L^(^n')\\p\\LP, 



which implies immediately 

\\p\\Lp{t) < C, 

with C independent of p, so our lemma follows. 

The key estimates on Vp and Vu will be given in the following lemma. 

Lemma 2.2 Under {\2.2^ . it holds that for any T <T*, 

sup {Vp\\l2 + llVwIlia) + f ||Vu||^idt < C. (2.5) 

0<t<T Jo 
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To prove Lemma 12. 2[ we need the following lemma (see [3] ) , which gives the 
estimate of V-u by divu and curlu. 

Lemma 2.3 Let u G H^{Q) be a vector-valued function satisfying u ■ n\Q^ = 0, 
where n is the unit outer normal of dQ. Then 

WuWh" < C'(||divtt||j:/s-i + ||curlu||j|/s-i + ||it||i^s-i), (2.6) 

for s > 1 and the constant C depends only on s and Q. 

Proof of Lemma \2.B . Multiplying p^^(/iAu + (/i + A)Vdivu — VP) on both 
sides of the momentum equations f lLip o. integrating the resulting equation over 
Q, one has after integration by parts 



dtj^r ' 2 



{(\\Yufdx+ / p ^{fiAu+ {fi + X)Vdwu-VPfdx 
Jn 



-p u ■ Vm ■ V X cnrludx + (2/i + A) / u ■ Vu ■ Vdivudx 
Jn Jn 



u-Vu- VPdx - Uf VPdx, (2.7) 
'n Jn 

due to Am = Vdivw — V x curlu. When u satisfies boundary condition (11. 4p or 

(II. Sp . we deduce from standard L^-theory of elliptic system that 



c||vp||i2 



f2.8) 



< CWfiAu + (/i + A)Vdivu||^2 - C\\VP\\l2 + C\\Vu\\l2 

< CWpAu + (/i + A)VdivM - VPllia + C\\Vu\\l2 

<C [ p'\pAu+ + \)VdiYu-VPfdx + C\\Vu\\l2, 
Jn 

due to p^^ > C^^ > 0. Lemma 12.31 yields that (12. Sp also holds for u satisfying 
boundary condition (11.60 due to the following simple fact by (11.60 : 

\\{2p + X)Vdwu\\l2 + ll/iV X curlM||^2 = \\pAu+ (p + A)VdivM||^2. 



Next, we shall treat each term on the righthand side of (12.70 under the boundary 
condition (II. 6p . since the estimate here is more subtle than the other cases, (11.40 
or (11.50 , due to the effect of boundary. 

Using (II. 6p and the facts that u x curlu = iV(| Vwp) — tt- Vw and V x (a x 6) = 
(6 ■ V)a — (a ■ V)b + (div6)a — (diva)6, one gets after integration by parts and 
direct computations that 



■ V)u ■ V X cnrludx 



/ curltt ■ V X ((tt ■ V)tt) dx 
Jn 

/ curlw ■ V X (m X cnT\u)dx 
Jn 



curl-updivudx — / curlu ■ T>{u) ■ curlwiix 



<c||Vw||i2(^)p(ii)IU. 



(2.9) 



and 



u ■ Vu ■ Vdivudx 



u^diu'njdiYudS 



Vu : Vu divudx H — / (divu) dx 
n 2 



< e\\V\\\l2 + C{e)\\Vu\\l2 {\\Vu\\l2 + \\V{u)\\lo^) 
due to the following simple fact: 



(2.10) 



dn 



u^diU^TijdwudS 



u^di{u ■ n)dwudS - 
an Jan 



u^u^ diUjdivudS 

an 

< C\\u\\l4i^Q^)\\dwu\\L2(an) 



u^u^ diUjdivudS 



< C\\Vu\\\2f^^)\\Vu\\Hii^) 



< C{e)\\Vu\\\.^^^^ + e\\VM\hi^n) + C{e) 



where fll.6p and the Poincare type inequality and the Ehrling inequality have 
been used. Similarly, 



u ■ Vu ■ VPdx 



an 



u^diU^UjPdS — / djU^diU^ Pdx — / u^didivuPdx 



an 



u'u^diUjPdS 



djU^diU^ Pdx — / (divw) Pdx — u ■ VPdwudx 
n Jn Jn 



< C||Vw||i2 + 



u ■ V Pdivudx 



< C||V'u|||2 + ||w|| re lldivwll r,3 II VP 



L6\\aiYU\\La II w -r\\L2 

5 1 



< C\\Vu\\l, + C\\Vu\\14V{u)\\1^\\Vp\\l2 

< ciivpiii^iiv^iiii^ + C\\Vu\\l. (PMIU^ + 1) + c, 

which yields also that 



(2.11) 



Ut ■ VPdx 



d 

dt Jn 
d 



Pdivudx 



Pfdivudx 



< 



. Pdivudx + u - V Pdivudx + (7 — 1) / P(divu) dx 
dt Jn Jn Jn 

^ [ Pdivudx + C\\Vu\\12{\\V{u)\\l^ + 1) 
dt Jq 

+C\\Vp\\UVu\\l2 + C. 



(2.12) 
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Substituting fl2:8|) - fl2J2|) into ([22]) gives that for e suitably small, 

^ j (^^\^u\^ + ^^{diYuf-P^lyv]dx + C4V\\\l, 
dt Jn \^ 2 / 

< C {\\Vp\\l + \\Vu\\l) {\\Vu\\l, + \\V{u)U^ + 1)+C. (2.13) 

It remains to bound the L^-norm of Vp. To this end, one can differentiate fll.ll) i 
and then multiply the resulting equation by 2Vp to get 

dt\Vp\^ + div(lVppw) + iVppdivu 
= -2(Vp)*VwVp - 2pVp ■ VdivM 

= -2(Vp)*P(m)Vp - 2pVp ■ VdivM. (2.14) 

Integrating (12.141) over Q yields 

dt (llVplliO < C\\Vp\\l4V{u)U^ + e\\V\\\l. + C(5)||Vp||i.. (2.15) 

Adding (12.151) to (12.131) . we deduce, after choosing e suitably small and using 
Gronwall's inequality, that (12. 5p holds. The proof of Lemma [2.21 is completed. 

Next step is to improve the regularity of p and u. We start with some bounds 
on derivatives of u based on above estimates. 

Lemma 2.4 Under the condition i\2.2^ . it holds that 

sup (||p^/V(t)||L2 + ||Vm||j^i) + / \\Vut\\l2dt<C, 0<T<T*. (2.16) 

0<t<T Jo 

Proof. Differentiating the momentum equations (II. ll) ^ with respect to t yields 

putt + pu-Vut — pAut — {p + XjVdiv Ut = —V Pt — PtUt — pufVu — ptU-Vu. (2.17) 

Taking the inner product of the above equation with ut in L^(f2) and integrating 
by parts, one gets 

\t+ [ P^tdx+ f {p\Vut\'^ + {X + p){divuty) dx 
2 dt Jn Jn 

= / Pfdivutdx — / p{ut ■ V u) ■ Utdx — / pu -V {\ut\^ + u -Vu ■ Uf^ dx 
Jn Jn Jn 

<C {\u\\Vp\ + \Vu\)\Vut\dx + C / {p\ut\^\Vu\ + p\u\\ut\\Vut\) dx 
Jn Jn 

+C / {\u\\ut\\Vu\'^ + \u\^\ut\\V^u\ + \u\^\\Vu\\Vut\)dx 
Jn 

= E^^- (2.18) 

i=l 

Noticing that by (12.51) and Sobolev's inequality, one has 

h < C(||u||loo||Vp||l2 + ||VM||i2) ||VMt||L2 

< C\\VUt\\L4^u\\Hi 

< 4'^ut\\l2 + C{e)\\Vu\\li. (2.19) 
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Similarly, 



h < C'||p^/^iit||L2 11^(11^6 II Vm||l3 + C||n||Loo||p^/2nt||i2||VMt||L2 

< C||p^/^Mt||L2||VMt||L2||VM||Hi 

< e\\Wut\\l2+C{e)\\p"^Ut\\l4yu\\''H., (2.20) 



and 



^3 < C'||M||i6||Mt||i:,6||VM||^3 + C||M^||L3||Mt||i6||V^w||L2 
+ C\\VUt\\LA\^u\\Le>\\u^\\L^ 

< C||Vut|U2 (||Vu|U2||Vm|U6 + ||V2m|U2) 

< £||V«,||i. + C(£)||Vti||^.. (2.21) 

We conclude from ^^-^M) that 

-— f pufdx + fi I \Wut\^dx 
2 dt Jq Jq 

< 6£||V«Ji. +C(5)||Vti||^. +C(.)||pV2„,||i.||V^i||l,i, 
which, together with Gronwall's inequality, implies that for e suitably small, 

sup ||p^/2ut(t) 11^2+ / \\Vut\\l2dt<C, 0<T<T*, (2.22) 

0<t<T Jo 

due to the fact that 

1 i i 

Po{x)'2ut{x,t = 0) = p^uo ■ Vuo{x) - PqQ e L^{^), 
which comes from the compatibility condition f ll.lSp . Moreover, since u satisfies 

pAu + (/i + A) Vdivw = put + pu ■ Vm + VP, 
fOD or (O) or dLEl) holds, 

similar to (12. 8p . one has 

\\V^u\\l2 < C{\\p^Ut\\L^ + \\u\\Lo.\\Vu\\L2 + \\VP\\L2) 

Hence, 

sup llVtillni < C. (2.23) 

0<T<T* 

Thus, Lemma 12.41 follows from f l2.22p and f l2.23p immediately. 

Finally, the following lemma gives bounds of the first derivatives of the density 
p and the second derivatives of the velocity u. 
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Lemma 2.5 Under the condition {\2.^ . it holds that for any q G (3, 6] 

sup WpWw^,, < C, 0<T <T*. (2.24) 

0<t<T 

Proof. In fact, fll.ll) i gives 

(iVpHt + div(|Vp|%) + iq- l)|Vp|MivM 

+ 9|Vp|^"=^(Vp)*P(M)(Vp) + gplVpl^-^p ■ VdivM = 0, 

which yields for the case that (11. 4p or (II. 5p holds, 

^llVplU, < C{\\V{u)\\l^ + l)||Vp|U, + C||VdivM||M, (2.25) 
and for the case that (II. 6p holds, 

|||Vp|U, < C(||I?(^)|U^ + IIVGIIm + l)||Vp|U„ (2.26) 

with G = (2/i + A)divM — P. Using the L^-estimate of elliptic system, we have for 
the case that (fLi|) or (fT3|) holds, 

IIV^mIIl, < C{\\put\\Li + \\u-Vu\\l, + \\VP\\l^) 

< C(||VMi|U2 + ||Vp|U, + l), (2.27) 

due to (I2.16p . When the boundary condition (II. 6p holds, noticing that (ll.6p yields 
that (V X curlw) ■ n = on the boundary dVL (see [2]), and (ll.ip i can be rewritten 
as 

VG = put + pu-Vu + pW X curlti, (2.28) 

we have 

VG ■ n\an = p{u ■ V)u ■ = -p{u ■ V)n ■ u\dQ,. (2.29) 
Therefore, (12.280 yields that G satisfies 

AG = div (p-Uf + pu ■ Vm) 
VG ■ n\9n = -p{u ■ V)n ■ u\dn. 

Using the L'^-estimate for Neumann problem to the elliptic equation, we have 

||VG|U, < G {\\put\\L^ + \\u ■ Vu\\l. + IIpImHIccq)) < C{\\Vut\\L^ + 1). (2.30) 

We deduce from fl2:25|) fresp. ^M)), (121271) fresp. UM)), fl216|) . and Gronwall's 
inequality that 

sup IIpIIvkm < G. (2.31) 

0<t<T 
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We complete the proof of Lemma 12.51 

The combination of Lemmas 12.41 and 12.51 is enough to extend the classical 
solutions of (p, m) beyond t > T*. In fact, in view of f l2.16p and f l2.24p . the 
functions {p,u)\t=T* = ^i^t^T*{p,u) satisfy the conditions imposed on the initial 
data (11.121) at the time t = T*. Furthermore, 

—fiAu — (// + A)V(divM) + \/P\t=T* = (P^t + ■ ^''^) — P^9\t=T*, 

whith g\t=T* G L'^{fl). Thus, {p,u)\t=T* satisfies (I1.13P also. Therefore, we can 
take {p,u)\t=T* as the initial data and apply the local existence theorem [HIS] to 
extend our local strong solution beyond T*. This contradicts the assumption on 



3 Generalization to the heat-conductive flows 

We can modify the previous proof to be fit for the heat-conductive flows. First, 
following the proof of Lemma 12.21 and noticing that 



Pt + u-VP + -fPdivu = (7 - l)KAe + (7 - 1)(|(Vm + Vu^f + X{diYuf), (3.1) 
one gets 



— f 7^|VmP + ^^-^^ — {divuy — Pdivudx 
at Jq 2 2 

+ / p'^{fiAu+ {p + X)Vdivu-VPfdx 



(3.2) 



<C{\\Vu\\lo. + 1)\\Vu\\12 + \ [ Ptdwudx\ 

Jn 

< Ci\\Vu\\L^ + l)||V«||i2 + e\\Vu\\jj^ + C(£)||V^||i2 
+ C(e)(l + \\Vu\\l^ + ||^||ioo)||VM||22, VO < £ < 1. 

Multiplying 9 on both sides of the energy equation, one has after integration by 
parts that 

dt I pe'^dx + 2K [ \Ve\'^dx <C\\9\\l^\\Vu\\l2. (3.3) 
Jn Jn 

Combining fl2.15p . (13.20 . (13. 3p and applying Gronwall's inequality, we easily de- 
duce that 

sup (||v..||i. + iivpiii. + Wphwi,) + r{\\vu\\l. + \\v9\\h)dt < a 

0<t<T Jo 

The higher regularity of (p, u, 9) can be obtained following the proof of [T8] . 
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